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ON THE IMPLICATIONS OF LATE INTERNAL DISSIPATION FOR SHALLOW-DECAY AFTERGLOW EMISSION AND 

ASSOCIATED HIGH-ENERGY GAMMA-RAY SIGNALS 

KOHTA MURASE 1,2 , KENJI TOMA 3 , RYO YAMAZAKI 5 AND PETER MESZAROS 3 ' 4 



ABSTRACT 

The origin of the shallow-decay emission during early X-ray afterglows has been an open issue since the 
launch of the Swift satellite. One of the appealing models is the late internal dissipation model, where X-ray 
emission during the shallow-decay phase is attributed to internal dissipation, analogous to the prompt gamma- 
ray emission. We discuss possible scenarios of the late prompt emission, such as late internal shocks, magnetic 
reconnection, and photospheric dissipation. We also consider the consequences of late dissipation and a two- 
component (early and late) jet model for the high-energy (GeV-TeV) emission. We study not only synchrotron 
self-Compton (SSC) emission from the early and late jets but also external inverse-Compton (EIC) emission, 
which is naturally predicted in the late dissipation model. For the latter, we perform numerical calculations 
taking into account the equal-arrival-time surface of EIC photons and show that the EIC component typically 
has a peak at ~ 1 — 100 GeV which may dominate over the SSC components. We demonstrate that very high 
energy gamma rays from both these components are detectable for nearby and/or energetic gamma-ray bursts, 
with current and future Cherenkov detectors such as MAGIC, VERITAS, CTA and HAWC, and possibly Fermi. 
Although the expected event rate would not be large, detections should be useful as a test of the model. Multi- 
wavelength observations using both the ground-based telescopes and the Swift and/or Fermi satellites are also 
important to constrain the models. 

Subject headings: gamma-rays burst: general — radiation mechanisms: non-thermal 



1. INTRODUCTION 

An understanding of the mechanism controlling the early 
X-ray afterglow emission from gamma-ray bursts (GRBs) 
has been one of the most debated issues since the launch 
of the Swift satellite. The canonical X-ray afterglow can 
be clas sified into three phases: th e st eep-decay phase (see , 
e.g.. | Kumar & Panaitescu 20001: [Yamazaki et al. 200*61 
Zhan g et al. 20061) . the shallow-decay (o r plateau) phase and 
the normal-decay phase (see, e. g.. iNousek et al. 2006 : 

O'Brien et al. 2006t iPanaitescu et al. 2006 : 

Zhang ^t al. 2006t IWillingale et al. 2007L and references 
therein). In particular, the shallow-decay phase is difficult 
to explain by the standard external forwa r d shock (FS ) 
model (see reviews, e.g., Meszaros 2006; Zhang 2007). 
Numerous models have been proposed so far to explain it. 
Most frequently discussed models are modified external 
shock models. One of the most popular interpretations 
involves either a continuous energy injection into the ex- 
ternal FS, where the long-lasting central engine energy 
output has a smooth decline oc T~ q , or short-lived central 
engine ejects shells with a steep power-law distribution of 
bulk Lorentz fact ors which can explain the shallow-decay 
emission (e.g.. iRees & Meszaros 19981 iDai & Lul998t 
IZhang & Meszaros 200 1 at IZhang et al. 20061)" Another 
version invokes a time-dependent microphysical scenario in 
the FS model, where changing parameters s uch as e e lead t o 
the observed shallow-decay emission (e.g.. Hokaetal. 20061) . 
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Dermer (2007) showed that the shallow decline may be 
explained if ultrahigh-energy cosmic rays are efficiently 
produced, by the recovery of an adiabatic relativistic blast 
wave after its radiative phase due to efficient photomeson 
losses and particle escape. The external reverse shock (RS) 
may also account for the shallow-decay emission, with an 
appropriate T distribution of the ejecta, if the RS emission 
dominates over the FS emission in the X-r ay band (e.g., 
iGenet et al. 2007t lUhm & Beloborodov 20071) . On the other 
hand, some authors have suggested explanations based on 
two- or multi-component jet scenarios. Panaitescu (2008a) 
showed that upscattering of FS photons by a relativistic shell 
can outshine the standard FS emission. Recently, Yamazaki 
(2009) proposed an alternative interpretation, where X-ray 
light curves are explained by the difference between the 
X-ray onset time and the burst trigger time. Among geo- 
metrical models, two co-aligned jets with different opening 
angles, i.e., wide and narro w jets, can also lead to the 
shallow-decay emission (e.g., Eich ler & Gran ot 2006), while 
a multiple-sub-j ets model was also proposed as one of the 
explanations dToma et al. 20061) . 

Another attractive interpretation is that X-ray emis- 
sion is attributed to lo ng-lasting internal dissipa- 
tion dGhisellini et al. 20071: iKumar et al. 2008bl) . This 
long-lasting dissipation or late prompt emission model 
can explain the chromatic behavior, which is not so easy 
to explain in modified external shock models where the 
optical flux would presumably track trends of the X-ray 
flux. This model is viable in the sense that the shallow- 
decay phase is not ubiquitous and some GRB afterglows 
are explained simply by the standard FS model when the 
late prompt emission is weak. Especially, some of the 
GRBs such as GRB 070110 have a plateau and a follow- 
ing steep decline, which strongly suggests that X-rays 
originate from the l ate internal dissip a tion rather than th e 
FS emission (e.g., iTroja et al. 20071 : iLiang et al. 2007bl) . 
This late internal dissipation model may also be consistent 
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with the existe nce of X-ray flar es in the early afterglow 
phase (e.g., iFalcone et al. 20071 IChincarini et al. 2010h . 
However, the situation is still inconclusive and unclear. 
For example, the lack of spectral evolution across the 
transition from the plateau to the normal-decay phase and 
the compliance of the closure relations in the normal- 
decay phase after the transition may rath er suggest FS 
models dLiang et al. 2007bHLiang et al. 2009b . 

On the other hand, recent novel results from the Fermi 
satellite have provided us with interesting clues for the 
mechanisms of GRBs. Especially, the onboard Large 
Area Telescope (LAT) has detected hig h-energy (> GeV) 
gamma rays from a fraction of GRBs ( Abdo et al. 2009a: 



lAbdo et al. 2009ft IXckermann et al. 20101) . Those detections 
have provided not only clues to the prompt emission mech- 
anism but have also provided the first detailed data about 
the high-energy afterg low emission, which had been expected 
for many years (e.g.. iDermer et al. 20001: [Sari & Esin 2001; 
Zhang & Meszaros 2001al) . In fact, late-time high-energy 
gamma-ray emission from GRBs such as 0809 16C, 090510, 
and 090902B has been attribut ed to afterglow emission rather 
than the prompt em i ssion dKumar & Barniol Duran 20101 : 
iGhisellini et al. 20101: iHe etal. 20101) . Various theoretical 
possibilities for the high-energy emission mechanism have 
also been discussed by numerous authors. The most 
widely discussed mechanisms are synchrotron and syn- 
chrotron self-Com pton (SSC) emission (see review, e.g., 
iFan & Piran 20081 and references therein). The exter- 
nal inverse-Compton (EIC) emission has been considered 
in some cases where seed photons come from flares or 
prompt emission ( e.g..lBeloborodov 20051 : [Wang et al . 20061 : 
iPanaitescu 2008b; M urase et al. 20101) . If protons and nuclei 
are accelerated up to very high energies, hadronic gamma-ray 
afterglows are also expected via th e photomeson production 
and ion synchrotron radiation (e.g., Bottcher & Dermer 1998; 
iPe'er & Waxman 20051: iMurase 2007l) ~ 

Despite recent progress in observations of high-energy 
gamma rays, the link between GRBs with shallow-decay 
emission and GRBs whose high-energy emission is detected 
by Fermi is uncertain, since due to the scarcity of simulta- 
neous detections with Swift it is unclear whether GRBs de- 
tected by ferm//LAT do have the shallow-decay phase or not. 
Also, the detectability with Fermi is limited at late times so 
that it is not easy to distinguish among the various models for 
shallow-decay emission. In this sense, Cherenkov detectors 
such as MAGIC and VERITAS may be more important. Al- 
though detections of > 10 GeV photons from distant GRBs 
become difficult because of the attenuation by the extragalac- 
tic background light (EBL), Cherenkov telescopes could pro- 
vide many more photons than Fermi when a nearby and/or en- 
ergetic burst occurs. Although very high-energy photons from 
GRBs have not been firmly detected so far, the future CTA 
and AGIS arrays would significantly increase the chances to 
observe high-energy GRB emission. 

Given a high enough detection rate of high-energy pho- 
tons by such observatories, high-energy gamma rays would 
provide very useful probes of the origin of shallow-decay 
emission. For example, Murase et al. (2010) demonstrated 
that EIC emissio n would be imp ortant to diagnose the prior 
emission model (Yamazaki 2009), which is one of the two- 
component (early and late) jet models. In this work, mo- 
tivated by the above prospects, we discuss theoretical pos- 
sibilities of late internal dissipation scenarios and investi- 
gate the associated high-energy emission. First, we review 



various late prompt emission scenarios, such as the late in- 
ternal shocks, dissipative photosphere, and magnetic dissi- 
pation scenarios. Second, we analytically study the high- 
energy gamma-ray emission expected in the internal dissi- 
pation model. In particular, we numerically calculate the 
EIC emission in detail, which plays an important role in 
two-component jet models such as the late prompt emission 
model. Throughout this work, cosmological parameters are 
set to Ho = 71 km s" 1 Mpc -1 , Qm = 0.3, and Q\ = 0.7, and we 
adopt the conventional notation Q = Q x x 10*. 

2. THEORETICAL POSSIBILITIES OF LATE PROMPT EMISSION 



Observationally, a good fraction of GRB afterglows 
show shallow-decay emission from T ~ 10 2 - 5 s to T ~ 
10 3 5 s (e.g.. lO'Brien et al. 200ft IWillingale et al. 20071: 
iLiang et al. 20 09), which can be expressed as 



L hP (T) oc 



(1) 



Here T a ~ 10 s is the break time when the shallow-decay 
phase ceases, a& ~ 0.0-0.5 and a st ~ 1.0-2.0. Many GRBs 
show the chromatic behavior, where optical and X-ray after- 
glows evolve in different ways, which tempts one to con- 
sider a two-component interpretation, i.e., X-ray and opti- 
cal emissions come from different emission regions. Some 
authors argued that the shallow-decay X-ray emission may 
be attributed to emission caused by internal dissipation sim- 
ilar to that of the prompt emission, while the normal-decay 
optica l emission is interpreted as an external FS compo- 
nent (IGhisellini et al. 20071: [Kumar et al. 2008b). For exam- 
ple, Ghisellini et al. (2009) successfully fitted X-ray and op- 
tical afterglows of various bursts in this picture. The isotropic 
radiation energy of late prompt emission, is typically 
££$ ~ (0.01-0.1) x 4 S ° B ~ 10 51 " 52 erg, where £™ B is the 
isotropic radiation energy of prompt emission. 

Such an interpretation seems strongly supported for a 
fraction of bursts. Some GRB afterglows show even a 
plateau rather than a shallow decay. For example, GRB 
070110 has the platea u, ctfl ~ 0.09, and the following sud- 
den decline, a st ~ 9 dTroja et al. 20071) . Such a behavior 
is very difficult to explain in the context of modified exter- 
nal shock models, which rather reflect variable activities of 
the long-lived central engine, although some of the issues 
can be sol ved if the emission is a nisotropic in the comov- 
ing frame (Bel oborodov et al. 201 lb . Then, the late prompt 
emission would presumably be attributed to long-lasting in- 
ternal dissipation. This interpretation seems consistent with 
the existence of flares, which are observed in about ~ 30 — 
50 % of GRBs ( IFalcone et al. 20071: IChincarini et al. 20101) . 
and flares may be attributed to accidental events of stronger 
intern al dissipation, in the late internal dissi pation sce- 
nario dGhiseflini et al. 20071: iKumar et al. 2008bl) . The dura- 
tion of flares ATf\ me is shorter than the observation time t and 
its flux enhancement is striking (the energy fluence is about 
10 % of prompt emission), which suggests that they origi- 
nate from temporarily strong late internal di ssipation by the 
long-lasting central engine (Iokaet al. 2005). One puzzling 
point is that the pulse width of flares increases linearly with 
time fchincarin i et al. 20101) . while those of prompt emission 
have various durations, showing no increasing pulse width 
during the burst. Possibly, this may reflect the behavior of the 
central engine. In the late internal dissipation model, one may 
expect some time variability in the apparently smooth X-ray 
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light curve of the shallow-decay emission, where it may also 
include some information on the central engine. However, 
at present, it is difficult to measure well even if the emission 
consists of numerous events of small internal dissipation (e.g., 
iGhisellini et al. 20071) . 

If the shallow-decay emission originates from long- 
lasting internal dissipation, what activity of the cen- 
tral engine could be responsible for it? For the 
central engine of GRBs, two possibilities have been 
most frequently discussed, accre tion of matter onto 
a black hole dKumar et al. 2008a|) or a fast rotating 
magnetar (e.g., lDai&Lu!998t IZhang & Meszaros 2001at 
Thompson et al. 2004t lYu et al. 2010h . In the latter scenario, 
the break time T a can be attributed to the spin down time. 



L hP (T) cx L P (T) cx 



const. 

j-2 



(T < T a ) 
(T a < T) 



(2) 



where Lp is the spin down luminosity. In this scenario, the 
outflow would initially be Poynting-dominated. 

In the former scenario, the shallow-decay behavior is at- 
tributed to the activity of the system of a black hole with an 
accretion disk, e.g., mass fall back accretion onto the cen- 
tral black hole. The break time t a can be interpreted as the 
end time of mass fall back. Kumar et al. (2008b) proposed 
that prompt emission is associated with the accretion of the 
innermost region of the progenitor star, whose angular veloc- 
ity is small, while the outer envelope with the larger angular 
velocity is responsible for the shallow-decay emission. The 
outflow luminosity is expected to be proportional to the mass 
accretion rate, and then the te mporal index ag is re lated to its 
behavior. After T a , we expect (Ku mar etal. 2 008a) 



L LP (r)«M B H(r)(x 



i+- 



T-T B 



2s — i r ac 



(3) 



where r ac is the accretion timescale and s ~ 0- 1. The above 
expression can explain both the rapid decline (when 7^ < T a ) 
and smooth transition (when T a < r ac ) at T ~ T a . 

There is also another interpretation of the behavior after T a 
and the origin of T a . If the late jet continuously decelerates, 
one expects a jet break in observe d light curves when Tlp 
becomes 0£p (IGhisellini et al. 2007b . 

Note that the late internal dissipation could also ex- 
plain the steep-decay emission just after the prompt emis- 
sion phase, although it is usuall y attributed to the high- 
latitude prompt emi ssion (e.g., iKumar & Panaitescu 20001: 
lYamazaki et al. 20061) . For example, in the collapsar sce- 
nario, Kumar et al. (2008b) suggested that the accretion of 
gas from the "transition" region between the core and the 
envelope, where the density has a steep decline, may lead 
to the steep-decay emission. However, such possibility that 
X-ray tails reflect the dying history of the central engine to- 
tally depends on models of the central engine, and the present 
situation is unclear since the apparent spectral evolution is 
also a ffected by the int rinsic spectrum of the prompt emis- 
sion dZhang et al. 20091) . 

At present, it is difficult to discriminate among these pos- 
sibilities from observations. Therefore, for the discussions 
below, we treat the temporal indices an and a st as just pa- 
rameters determined from observations. Also, we assume that 
the long-lasting internal dissipation occurs according to Equa- 
tion (1) without specifying the central engine. 

2.1. Late Internal Shock Scenario 



In the classical scenario, the prompt emission is ex- 
plained by electromagnetic radiation from electrons accel- 
erated at internal s hocks that occur in the optically thin 
relativistic outflow (Rees & Meszaros 1994). Flares may 
also be explained similarly, where X-ray and/or ultravio- 
let photons are produced by relativistic elect rons acceler- 
ated at late internal shocks dFan & Wei 2005b . The bulk 
Lorentz factor responsible for flares is ofte n thought to b e 
smaller than that of prompt emission (e.g., iJin et al. 2010), 
and then, applying this scenario to late prompt emis- 
sion, the typical collision radius is estimated as r,- s» 
2r2 p <57\ ai - ~ 1.5 x 10 15 cm (r LP /5) 2 (5r var , 3 (l +z)~ l . (How- 
ever, i n some models such as the fast rotating magnetar 
model dThompson et al. 20 04). the bulk Lorentz factor of the 
late jet may be much larger.) 

Let us consider the two-shell collision between fast and 
slow shells with Tf and T,, respectively. The relative Lorentz 
factor between the shells is r LP . sh « (T f/T s + T s /T f )/2 — 5 
and the average Lorentz factor of random internal motions 
is r LP , ls a ( V /T77r;+ y/TJT})/2 ~ 0.5 for Tf ~ .13 and 
r s ~ 2 (where the Lorentz factor of the merged shell is Tlp ~ 
y/T fT s ~ 5). Then, the electron injection Lorentz factor is 
obtained as 



7e, 



(ii 



1)- 



0.92 x 10 



LP.is " 



1 



0.5 



(4) 

where e e is the fraction of the internal energy transferred 
to non-thermal electrons and f e is a number fraction of 
accelerated electrons. Introducing which is the frac- 
tion of the internal energy transferred to the magnetic field, 
the comoving magnetic field is estimated as B ~ 3.5 x 



10 3 G 



(r L p.i,+3/4xr L 



-i) 



(9/8) 



1/2 



4^i L M9 r j5( r Lp/5) l > and then 



the observed synchrotron peak energy is 

"(r LP , is -i) 5 / 2 (r LP , is +3/4) 1/2 



;0.17keV 



(3/8V2) 



x fe^l e l-l e B,-l^k'.49 r i.l5(^ + z ) 



1/2 ,1/2-1 



(5) 



The observed X-ray emission shows a hard spectrum in the X- 
ray band, F^p cx E~ l , which is attributed to synchrotron emis- 
sion by relativistic electrons with the spectral index of p ~ 2. 
Accelerated electrons are expected to be in the fast cooling 
regime. 

In the case of prompt emission, the classical scenario has 
several problems in explaining observations. One of them 
is that, if the radiation mechanism is synchrotron emission, 
explaining the low- energy spectral inde x of /?; ~ 1 is not 
easy (but see, e.g ., iDerishev et al. 200lt iNakar et al. 20091 : 
Daign eetal. 20111) . But, since the spectrum in the far- 
ultraviolet range is not observed, it is uncertain whether the 
late internal shock scenario suffers from the same issue. An- 
other one is the efficiency problem that energy transferred to 
prompt emis sion often seems la rger than the afterglow kinetic 
energy (e.g., [Zhang et al. 20071) . The efficiency seems worse 
due to weaker collisions in the late internal shock model, but 
it is observationally unclear whether the shallow-decay emis- 
sion has the same problem. 

2.2. Dissipative Photosphere Scenario 

In the previous subsection, we discussed the late internal 
shock model based on the analogy to the prompt emission. 
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However, the prompt emission mechanism itself is still 
under debate, and many possibilities have been suggested. 
Another popular scenario of the prompt emission is the 
photospheric emission model, where quasi-thermal emission 
comes from around the pho tospher e (tj = n e oj(rjlY~) 



1) (e.g., [T hompson 



Meszaros et al. 200 



T 



1994; 



Meszaros & Rees 2000; 



Rees & Meszaros 2005; 



Pe'er et al. 20061). Although there are various ver - 

sions (e.g.. lIokaetal. 20071 lBeloborod ov20T(it EokaJOlO), 
we h ere consider the dis si pative photosph ere sce- 
nario dRees & Meszaros 20051 iPe'er et al. 20061) . where 
internal dissipation (via e.g., internal shocks or magnetic 
reconnection) occurs around the photosphere. For the late 
jet making late prompt emission, the photospheric radius is 
written as 



r p h : 



e r 47rr?p»ipC 3 



;9.4 x 10 13 cm (C/10e f )L L p,48(r LP /5)" 



(6) 



where e r is the ratio of the radiation energy to the kinetic en- 
ergy carried by cold baryons, and £ e is the ratio of the number 
of electrons to the number of baryons, taking into account 
the possibility of copious pair production via internal dissipa- 
tion. The photospheric radius thus obtained would generally 
be above the typical radius of a Wolf-Rayet star. Also, es- 
pecially when the late jet is baryon-rich compared to that for 
prompt emission, the photospheric radius is likely to be lo- 
cated above the coasting radius. The comoving temperature 
at the photospheric radius is 



kT ~ 30 eV 4?48'-ph 1 / i3.5(r L p/5)- 1/2 , 



and the observed typical energy is 



7 b ~ 48 keV f 5/12 r -2 ' 3 



10 



1/6 



(7) 



(8) 



The observed X-ray emission shows a hard spectrum in the 
X-ray band, Flp oc E" 1 , which requires some process such 
as Comptonization by nonthermal electrons produced via in- 
ternal dissipation around the photosphere. In this scenario, 
the variability timescale would be relatively short, 5r v ar ~ 

21 s r,.i3 5(Tlp/5) 2 (1 +z), although dissipation itself may last 
for a longer time. 

As mentioned before, flares would also be caused by ac- 
tivities of the long-lasting central engine. However, it might 
not be easy to explain flares in this scenario. Flares seem to be 
caused by occasional larger dissipation of relativistic outflows 
but those with Lorentz factors whose values are larger than the 
values of prompt emission but may be sm aller than the values 
of late prompt emission, rg are ~ 10-50 dJin et al. 20101) . On 
the other hand, light curves of flares often show the expo- 
nential decay after the peak, which suggests relatively large 
emission radii of ~ 3 x 10 14 cm rg are j Arfl are ,2(l +z)~ l , if 
the decay of pulses is attributed to high-latitude emission. The 
typical radii seem above the photospheric radius, unless the jet 
is largely pair dominated. Here, one should keep in mind that 
the photospheric scenario and other scenarios are not mutually 
exclusive. For example, internal shocks may occur well above 
the p hotospheric radius as we ll as around the photospheric ra- 
dius dRees & Meszar os 2005b . Magnetic reconnection is also 
one of the possibilities. 



2.3. Magnetic Dissipation Scenario 

Relativistic jets launched by the central engine may be ini- 
tially Poynting-dominated. If the outflow is still Poynting- 
dominated at the emission radii, without significant conver- 
sion into the kinetic energy, magnetic dissipation rather than 
shock dissipation of the bulk kinetic energy may lead to pro- 
duction of nonthermal particles. Although detailed scenar- 
ios for this are still unavailable due to lack of our knowledge 
on mechanisms of magnetic dissipation and associated parti- 
cle acceleration, prompt and/or late prompt emission may be 
produced by internal dissipation of a significant fraction of 
the magnetic energy in the outflow. For example, Lyutikov 
(2006) argued that magnetic dissipation may occur around 
the radius where the MHD approximation breaks down, if the 
outflow is extremely magnetized. On the other hand, mag- 
netic fields are distorted by internal shocks at r,- ~ 10 13-15 cm, 
whic h may eventually lead to efficient magnetic reco nnec- 
tion dZhang & Yan 201 UlMcKinnev & Uzdenskv 20101) . 

In this work, just for demonstrative pu rposes, we apply the 
jets-in-a-jet model (Giannios et al. 2009| for the late prompt 
emission. In this scenario, the magnetic reconnection in a 
jet leads to many mini-blobs with relative Lorentz factors of 
~ yfa. Indeed, radiation from such mini-blobs can reproduce 
highly variable light curves of promp t emission, though there 
remain potential problems (see, e.g., lLazar et al. 20091) . As- 
suming that the late jet has r LP ~ 5 with the magnetization 
parameter of a ~ 30, the dissipation radius is estimated as 
n ~ 3.8 x 10 15 cm (r LP /5) 2 a 1 . 4 <5r var , 2 (l+z)- 1 . 

The typical electron Lorentz factor is estimated as 



e e Vv— ^ 0.92 x lOVj 72 e. 



(9) 



The magnetic field in the downstream blob can be fi ~ 
5.2 x 10 2 G L}Jl s rj\ 5 5 (T LP /5)~ l (note that it does not have 
to be Poynting dominated since a significant fraction of the 
magnetic energy is dissipated there). Then, the typical syn- 
chrotron peak energy is estimated as 6 



.0.13keV £ 2 _ lt r 3/ 4 2 42 8 r;, 15 . 5 



(10) 



which seems consistent with observations. A hard spectrum, 
Flp oc £■"' , may be attributed to synchrotron emission from 
nonthermal electrons accelerated at shocks caused by mini- 
blobs. Note that electrons are typically in the fast cooling 
regime, since 7<,, c ~ 0.14-Lg 1 4S ri,\5.5(Ti J p/5) i (of course, the ac- 
tual electron Lorentz factor should be larger than unity). Here, 
the lower bulk Lorentz factor with the lower magnetization is 
assumed for the late jet compared to the case of prompt emis- 
sion, but it may not be the ca se in some models s uch as the 
fast rotating magnetar model (Metzger et al. 2010h . In order 
to have the synchrotron peak of ~ 0. 1 keV, not only large radii 
but also low values of e e might be required in this scenario. 
For example, if a ~ 10 5 and r, ~ 10 17 cm, e e ~ 10" 3 and the 
short variability time are expected, depending on the scenario. 



3. ASSOCIATED HIGH-ENERGY EMISSION 

Next, we consider consequences of the late internal dis- 
sipation model for high-energy emission. One can consider 
two possibilities of high-energy emission. One is EIC emis- 
sion produced by electrons accelerated at the external shock 

6 If one applies this model to GRB prompt emission, we have E b ~ 
0.57 MeV ^.^Ous.sO+zr'- 
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caused by the early jet, which is responsible for the prompt 
emission and the observed standard afterglow component. 
Late prompt photons from inner radii are naturally upscat- 
tered in the late internal dissipation model, and predictions 
are not sensitive to details of late internal dissipation models. 
In this paper, we especially discuss this possibility in detail 
(see the next section). The other is the high-energy emission 
from the emission radius at which internal dissipation occurs, 
e.g., SSC emission from the late jet. Obviously, predictions 
of high-energy emission depend on each scenario. This pos- 
sibility is also discussed in this section. 

3.1. High-Energy Afterglow Emission 

First, we discuss EIC emission caused by interactions with 
late prompt photons and electrons accelerated at the external 
shock of the early jet producing prompt emission. As demon- 
strated below, this EIC emission is useful as a test of the late 
internal dissipation model. We here give analytical consider- 
ations, but more detailed results with numerical calculations 
are provided in the next section. 

We can think that one of the two components is the stan- 
dard afterglow component from the early jet. For an adi- 
abatic relativistic b last wave expanding int o the interstel- 
lar medium (ISM) (Blan dford & Mckee 1976b . we obtain the 
bulk Lorentz factor as 



7ej£ ~ TT{lehe.m) for J e , m < -f e < -f e , c and J e j£ 

TT{le,c/le.mY P+l {le/l e ,cY P for J e > Je.c, the resulting EIC 
spectrum in the Thomson limit is expressed as 

( E 2 ~P> (E<E" IC ) 
EF EIC (E) cc I E^/ 2 (Eg IC <E< E^ IC ) (15) 

[ ecw (E c mc < E) 

where tj ~ (ojN e i 'A-ttR 2 ) is the Thomson optical depth, M e 
is the number of electrons, /3i is the low-energy photon index 
of late prompt emission, and 



(ID 



T(T) ~ 44 £l! 5 ]n- U % y \l + z) 3/ \ 

and the external shock radius is estimated as 

R(T) ~ 2.3 x 10 17 cm S^ 1,A T^\l + zT l/ \ (12) 

where is the isotropic kinetic energy of the ejecta and n is 
the ISM density. 

Electrons would be accelerated at the external FS. The in- 
jection Lorentz factor of electrons is estimated as 

lem ~ 2.3 x \Q\ ef ,- 2 f;}{g p /g 2A )£ l k ^n;; l *T- y \\+zf\ 

(13) 

where g p = (p- l)/(p — 2) and p is the spectral index of 
FS electrons. Here e e f is the fraction of the internal en- 
ergy of the shocked ISM transferred to non-thermal elec- 
trons at the external FS, and f e f is the number fraction of 
elec trons injected to the acce leration process at the external 
FS (lEichler & Waxman 20 05). The cooling Lorentz factor of 
electrons is estimated by fd yn = fcooi> an d we nave 



7e 



■ 6.2 x 10 4 e 



dVVd+^d+y)" 1 d4) 



B/,-3°jfc,53 "0 



where fd yn = A/c ~ (4/k)TcT is the dynamical timescale, 
f coo i is the electron cooling timescale, and Y is the to- 
tal C ompton Y parameter. He re, k is set to 4 in this 
work (Panaitesc u & Kumar 20041) . and eg/ is the fraction of 
the internal energy of the shocked ISM transferred to the 
downstream magnetic field. In the slow cooling case (7^ < 
7e jC ) with a constant Y, the steady electron distribution is 
dJ\f e /d-f e oc -f~ p for 7 tv „ < j e < 7^. and dN e /d^ e oc % p ~ 1 for 
7e,c < 7e- In the fast cooling case (7^ < j e , m ) with a con- 
stant Y, the steady electron distribution is dJ\f e /dy e oc 7J 2 for 
7e,c <le< le,m and dM e / 'dj e oc 77" 1 for 7^ < 7 e . 

These electrons upscatter late prompt photons at the vicin- 
ity of the FS. The expected EIC luminosity is very roughl y 
written as L E ic ~ min(y E icL LP ,L e ) (e-g-, iFan et al. 2 008). 
where J^ic is introduced as the ratio of the EIC energy 
flux to the seed photon energy flux. In the slow cool- 
ing case, noting that F E ic(£) ~ J d^ e ^F LP (j e ,E), where 



"EIC ~ le.m^ 



; 5.2MeV£ 6 ikev4-2/«/l — I &«"■ 



J EIC ' 



g2A 



.95 GeV E D A ^/.-^V' 4 [ 



1/4 -1/4 



■"t.53"0 



1/4 



T 3 

l+z 



l+Y 



-3/4 

(16) 



(17) 



The contribution below E™ sc mainly comes from interactions 
between electrons with ~ j e>m and photons with E < E , while 
the contribution in the range E™ sc < E < E$ sc comes from 
interactions between electrons with j e m < -f e < j e c and pho- 
tons with ~ E b . The EIC flux at E EIC is also estimated from 
E EIC F£ IC sa xYeic(E c F c ), where the EIC Compton Y parame- 
ter, Yeic. is introduced as the ratio of the EIC energy loss rate 
to the synchrotron energy loss rate (and it is different from 
J^eic)- Here, x < 1 is a factor coming from the fact that the 
EIC emission is anisotropic. 

The Thomson limit has been consider so far. However, the 
Klein-Nishina (KN) effect o ften becomes importan t at suf- 
ficiently high energies (e.g., iGuetta & Granot 20031) . In our 
case, there are three characteristic energies. 



E"^^T 7ejn m e c 2 /(l+z) 
E l KN ~T"/ ex m e c 2 /(l+z) 
Ei N ^2T 2 m 2 c 4 /E b /(l+z) 2 



(18) 
(19) 
(20) 



When the KN effect becomes important, the EIC spectrum has 
breaks. Here, let us introduce £kn.i as the first break energy 
due to the KN effect. When the KN break exists above the 
EIC peak, instead of Equation (15), we have 



EF E ic(E) oc 



E 2 ~P< (E < £™ c ) 

E (3- P )/2 {E m ic <E<E£ lc ) 

E (2- P )/2 (££ IC <E< £kn,i) 

EPr-P (£kn,i < E) 



(21) 



Here 



- F 



: 9.9 TeV « ^'ft^d ^ . 

(22) 

This case is typical for our adopted parameters, and the EIC 
emission at E > £kn.i is dominated by radiation from elec- 
trons with 7 e ~£/rm f c 2 (l +z) interacting with seed photons 
with the energy of ~ Y 2 m 2 e c A /E{\ +z) 2 via the Thomson scat- 
tering. 

If the KN break appears below E^q, we obtain 



E 2 ~P> (E < £* c ) 
EF E1C (E) cx { E^l 2 {E% lc <E< E KNA ) 
EP<-p (E kna <E) 



(23) 
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where 

£KN,i=£™^690GeV e^_ 3 ^ 



-1/4-3/4^-1/4 



*.53 "0 



'3 



(1+z) 



-3/4 _ 



1+F 
(24) 



If 7 ejBI and/or are too large, one expects the deep KN 
regime. In this case, we have 

C £' 2 "' 9 ' (E < £ K n,i) 

£F EIC (£) a 1 ePi-p +1 (Ekna <E< £kn. 2 ) (25) 

[ E*-p (£ KN .2 < E) 

Here, 



: 5.1 GeV e e/ ,_2# fe w 



^.4)^ 4 3 «o 1/4 ^ 3/4 (i^r 1/4 

(26) 



and £kn,2 = r7 CiC ra ( ,c 2 /(l +^) 2 is the second KN break. This 
EIC spectrum is anticipa ted in the prior em ission model for 
shallow-decay emission (Mura se~et al. 20101) . 

We are interested especially in cases where the EIC flux 
exceeds the afterglow SSC flux. For this purpose, we next es- 
timate the SSC flux. The SSC e mission has been st udied by 
many authors (see reviews, e.g., iFan & Piran 20081) . so that 
we here discuss it just briefly. The char acteristic energies 
of the SSC emission are o btained as (e.g.. ISari & Esin 2001t 
IZhang & Meszaros 2001bT) 



,1/2 
ef c e/,-2 e fl/,-3 



£ s ffl sc -4.7 keV {gp/giAff^e'e 



££ sc ~490 GeV 



(27) 



l+Y 



-7/2 
£ Bf-3 



x^W /4 d^r 3/4 - (28) 

For the slow cooling case that we are interested in, the SSC 
spectrum in the Thomson limit is expressed as 

f E^ (E<E£ SC ) 
EF S sc(E) oc I £C3-p)/2 (£« c < e < E c ssc ) (29) 

[ £(2"P)/2 ( £ c sc < £) 

Note that only the first SSC component is important, since the 
second SSC component is typically negligible due to the KN 
suppression. The energy flux at the SSC peak is also evaluated 

as 

p-3 



E ssc F ssc - 1-2 X 10" 7 GeV cm" 2 s" 1 F ssc 



l+Y 



l L,21.S 



_2J_ | f 2-p p-l P-2 C2 h 
I Jef t ef.-2 t Bf,-3 c k,53 n O 



gp 
82A 



by which we can normalize the SSC spectrum. As a re- 
sult, the EIC flux and SSC flux are roughly related as 
^eic^eic/^SSC^SSC ~ ^eic/Tssc- See Section 4 for detailed 
discussions on the relative importance of each component. 

The KN effect may become important in the cases we con- 
sider here. For our typical parameters, the KN break is located 
above EZ SC , which is 



- 2 mlc A 



(1+z) 2 £ c 

lAToYel^S^nl^il+zT 



l+Y 



(31) 



In general, SSC spectra can be complicated and co nsist of sev- 
eral breaks (Na kar et al. 20091: IWang et al. 20ldb . Hence, we 



numerically calculate the SSC emission taking into account 
the KN effect. 

We can also calculate light curves, once the dynamical evo- 
lution of the blast wave is given. In the next section, we show 
the resulting light curves of the EIC and SSC emission. Note 
that the temporal behavior would change after the jet break 
time of Tj - 10 5 s dRhoads 1999t ISarietal. 19991) (where 
T9j ~ 1), but throughout this work we focus on the behav- 
ior before tj. 

3.2. High-Energy Late Prompt Emission 

High-energy emission is expected from the late jet itself, as 
mentioned before. For example, one can expect SSC emis- 
sion as well as synchrotron emission if electrons are accel- 
erated in the magnetized region. If protons are also acceler- 
ated up to very high energies, hadronic gamma rays are pro- 
duced via photomeson and photopair production, and proton- 
synchrotron radiation. Predictions depend on models, which 
are quite uncertain. Hence, in this subsection, we just provide 
analytical considerations on several interesting cases. 

In the photospheric scenario, the injection Lorentz fac- 
tor of electrons should be 7 e m ~ 1 to produce the hard 
spectral component of Flp oc E~ l by Comptonization (e.g., 
iThompson 19941: lloka et al. 200% . The IC spectrum may be 
extended up to high energies, but high-energy photons cannot 
avoid attenuation by pair-production. Murase & Ioka (2008) 
showed that the pair-production break (or cutoff) should be 
around £ cut w ^■T]jpm e c 1 /(1+z) in the pair-photospheric sce- 
nario (for fit, = 2), which suggests that high-energy emission 
above GeV is not expected in the one-zone case. In the multi- 
zone case, relativistic electrons may be produced at outer 
radii. For example, internal dissipation may occur above the 
photospheric radius, leading to EIC emission with the typi- 
cal energy of ~ r i l em E b ~ GeV. However, we will not discuss 
here such more complicated possibilities. 

In the magnetic dissipation or late-internal shock scenarios, 
the typical emission radii are much larger, so that it is easier to 
expect high-energy gamma rays that escape from the source. 
Here, as a demonstrative example, we consider the SSC emis- 
sion in the jets-in-a-jet model described in the previous sec- 
tion, which is sufficient for our purposes in this work. First, 
in the Thomson limit, the typical SSC energy is estimated as 



■ 210 MeV ei^L^lr-'il +zf 



(32) 



Introducing the Compton Y parameter, Ylp, the SSC flux at 



"SSC 



is written as 



J SSC r SSC ~-*LP 



i \ ph ti P 



-2.5 x lO" 6 GeV cmV 1 YufL ^ \ +Z .(33) 

In the Thomson limit, Ylp can be approximated as Tlp ~ 
i+^/i+4e c /eB ^ ^ e ff ect ma y actually become impor- 
tant at sufficiently high-energies. When the KN break exists 
above E$ sc , we have 

< 2 "A (E<E b sc ) 
EF SSC (E) cx < E 2 ~P" (E b sc <E< #kn,i) (34) 

The KN break is given by £ K n.i = (T em m e c 2 /(l + z)E b ) 2 E b ~ 
1.3 TeV (r L p/5) 2 e- 2 1 a- 1 4 /2 L-| 4 / 8 2 r ; - 15 . 5 (l+z)- 1 . If £ K n,i < 
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£g SC , the spectrum is in the deep KN regime, and we obtain 



EF SSC (E) oc 



£ 2 "A (E < £kn,i) 

r e m7e I m« e C 2 /(l + Z) 



(35) 



where £kn.i = 

12 GeV e e -i(rLp/5)f7i.4(l+z) _1 - In the fast cooling 
case, the resulting spectra can be more complicated espe- 
cially when Flp in the Thomson limit is so large that the 
distrib ution of electrons is affected by the KN effect (see, 
e.g., iNakar et al. 20091; iBosniak et al. 20091) . But the above 
expressions are reasonable for moderately small values of 
*lp. 

In the late internal dissipation scenario, the pair- 
creation process is crucial for high-energy gamma-ray emis- 
sion. The optical depth for pai r production is esti- 
mated as (e.g., iLithwick & Sari 200 it iMurase & Ioka 20081; 
iGupta & Zhang 20081) 

0-i 



'77 



0.l(7 T l 



L b 



/ (l+zTEE b 
47rrfT em c(l+z)E b { T 7 - m my 



(36) 



where I is the comoving width. Assuming / ~ r,-/T em , 
the pair-production break (or cutoff) is estimated as E cut ~ 

I 1 20 P 2/3 3-2B 



8.5 GeV (L» p48 ) ^ 5 (r LP /5)^a^ (E b A keV )^(l +z)^ 
in the magnetic dissipation model. At energies higher than 
this energy, the spectrum is suppressed or may have a cutoff. 

In Figures 1 and 2, we show SSC spectra which are calcu- 
lated analytically using Equation (1341 . The parameters are 
described in the caption of Figure 1. The p air-productio n 
opacity is taken into account by 1/(1 + r 77 ) (Baring 2006). 
From Figures 1 and 2, we see that it is difficult to detect high- 
energy late prompt emission from distant bursts, but is possi- 
ble for nearby bursts. The SSC peak is expected around the 
GeV range, which may be reached by Fermi if GRBs occur at 
z < 0.7. Very high energy gamma rays above ~ 30 GeV are 
also detectable with the future CTA, although its detectablity 
depends on the pair-creation opacity both inside and outside 
the source. In Figures 1 and 2, the KN break is seen around 
TeV but the attenuation by pair-creation masks it. The light 
curves of the high-energy gamma rays basically follow the 
observed X-ray light curve of the late prompt emission. 

One can also calculate the SSC emission in the late in- 
ternal shock scenario similarly to how it was done in the 
previous paragraph, by changing parameters. In this para- 
graph, we briefly discuss the hadronic emission, although 
detailed studies are beyond the scope of this work. In the 
late internal shock scenario, not only electron but also pro- 
tons may be accelerated up to very high energies. Even in 
the magnetic dissipa tion scenario, protons may be acceler- 
ated dGiannios 2 010). although a large baryon loading may 
not be expected. Hadronic emission in the late internal dissi- 
pation model was considered and discussed in Murase (2007). 
The Lorentz factor of the late jet might be relatively small, 
and the late jet might be more baryon-rich compared to the 
early jet making prompt emission. Then, as shown in Murase 
(2007) and Murase & Nagataki (2006), copious soft photon 
fields in the late jet lead to a high meson production efficiency 
given by 



fi 



in 



1.4- 



L LP/ 



<E P /E b /-\ 



(37) 



1 keV 

where the multi-pion production effect (which is a factor of 
three) is taken into account. The expected neutrino flux is 





-6 - 


'{/) 




CM 




cm 


-7 - 


> 




<D 




CD 




u_ 


-8 7 


LJJ 




CD 
O 


-9 - 



-10 



SSC 

Syn 

Fermi/LAT 

CTA 




"log(E y [G 2 eV]) 



FIG. 1. — Spectra of synchrotron and SSC emission in the magnetic dis- 
sipation model for late prompt emission at T = 10 3 s. The source redshift 
is taken to z = 1. Assumed parameters are: Lb\t> = L e \fi = 10 48 ' 5 erg s , 

a = 10 1 ' 4 , IYp = 5, and r, = 10 15 ' 75 cm. The thick solid curve represents an 
SSC spectrum taking into account attenuation by pair-creation both inside 
and outside the source. An SSC spectrum shown as the thin sold curve in- 
cludes only the source attenuation, while the thin dotted curve spectrum does 
not include either of them. The Fermi/LAT and CTA sensitivities (with the 
duty factor of 30 %) are also overlayed (CTA Consortium 2010). The LAT 
sensitivity curves in the sky survey mode are used for the long time obser- 
vations, although the possible continuous observations by LAT may improve 
the detectability by a factor of 3-5 (e.g., Gou & Meszaros 2007). 
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FIG. 2. — Same as Figure 1, but T = 10 s and ; 



= 0.3. 



comparable to or maybe larger than that of prompt emis- 
sion, since the meson production efficiency is high while 
the total radiation energy of flares or l ate prompt emission 
is ~ 10 % of that o f prompt emission (Falcone et al. 20071 
IChincarini et al. 20101) (and the kinetic or magnetic energy 
may be larger if the radiation efficiency is low). Hadronic 
gamma rays are also expected as well as neutrinos. The ef- 
ficient pair-production in the source induces electromagnetic 
cascades. Assuming an E~ 2 spectrum for cascade emission, 
the gamma-ray flux is crudely estimated as 

1+z 5L CR 



EF n , 



1 .6 x 10" 6 GeV cnrV 1 L C r 49 - 2 ° ' " 



(38) 



TZ d 2 

/V « L;2 7.5 

where TZ is the conversion factor from the total energy amount 
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of protons into the energy amount of protons per energy 
decade, 1Z ~ 20 for p=2. Although the detailed calculation 
is beyond the scope of this work, this suggests the potential 
importance of hadronic gamma rays for nearby GRBs. 

Although predictions for the high-energy late prompt emis- 
sion are model-dependent, once high-energy gamma rays are 
detected in sufficient amounts, they would be useful to dis- 
tinguish between the various uncertain mechanisms of late 
prompt emission. 

4. NUMERICAL RESULTS OF HIGH-ENERGY AFTERGLOW 
EMISSION 

In the previous section, we gave analytical estimates of EIC 
and SSC emission produced by relativistic electrons acceler- 
ated at the external FS. As we discussed, the EIC emission 
does not depend on details of late internal dissipation mech- 
anisms, and is useful as a good probe of the different scenar- 
ios. In this section we calculate numerically the EIC emis- 
sion, which provides significantly more accurate results than 
the analytical estimates. This is because: (1) the EIC emis- 
sion is anisotropic, which leads to suppression by a factor of 
x; (2) the KN suppression becomes important at high energies 
above £kn,i; (3) the influence on the electron distribution is 
complicated if the EIC/SSC cooling is efficient and the KN 
effect is relevant. 

In order to calculate the EIC emission, we need to consider 
the equal-arrival-time surface of upscattered photons. The ex- 
pression for the EIC emission is written as (see Appendix A) 



Fe\c{T) = -ctt 



\-2y+2y 2 + 



2(1-0 



(i+r 2 2 y 



where y ■ 



and ^ir 1 ^ The scatter- 

^ 7.1 -v..m..r- 



2(l-cos9)7 r e(l-0 "» — IT-i.m^c 1 

ing angles 9 and 9 of EIC photons are measured in the central 
engine frame and the comoving frame, respectively. The func- 
tion G(e) represents the spectral shape of seed photons with 
energy e in the comoving frame (e.g., e = (1 +z)E h /2T). In 
the case of a broken power-law spectrum for late prompt emis- 
sion, it is G(e) = (e/e h )~ Pl+l for s < e h and G(e) = {e/e h ) Ph+l 
for e < e, respectively. 

The input parameters required for the calculations are basi- 
cally determined by afterglow observations at X-ray and op- 
tical bands. We set typical parameters following Ghisellini 
et al. (2009). As for the electron distribution, we exploit 
the standard ext ernal FS model (e.g., Meszaios & Rees 19971: 
ISari et al~9 98) and adopt the following fiducial parameter 
set: S k = 10 53 - 5 erg, n = 1 cm" 3 , e ef = 10" 2 , e Bf = 10" 3 and 
p = 2.4. 

For the late prompt emission, the seed photon spectrum is 
assumed to be a broken power-law spectrum with /3/ = 1 for 
E < E b and (3 h = 2.2 iorE b < E, with E' b = 10 2 - 10 2 ' 5 eV. The 
break time is set to T' a = 10 3 s and the late prompt luminosity 
at T' a is taken as L b LP \ T > = 10 48 - 10 48 5 erg s" 1 (which means 
£hP,x/tef£k ~ 1)- The temporal indices are set to an = 0.2 
and a st = 1.5. Also, assuming r,- = 10 13 ' 5 cm, the high- 
energy cutoff due to pair creation is determined from E cnt ~ 

I _!_ 20 , 20 3 _2/3 

8.5 GeV (Lf P4S )"^r . ^ ,(r em /25)^(£ fc , k v )^(l +z)?r (e.g.. 




1 2 

log(E y [GeV]) 

FIG. 3. — Gamma-ray spectra of EIC emission in the late internal dis- 
sipation model for GRB afterglows, caused by Compton scatterings of X- 
ray photons by electrons accelerated at the ext erna l shock. Calculation are 
numerically performed according to Equation )39h taking into account the 
equal-arrival-time surface. The observation time is set to T = 10 3 s and the 
source redshift is taken as z = 0.3. Relevant parameters for the late jet are 
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iGupta & Zhang 20081: iMurase & Io ka 2008) with the atten- 
uation factor of I/IT + T-y^(E)] dBaring 2 006). and the 



erg s _1 , r a ' = 10 3 s, E" = 0.1 keV, a„ = 0.2, and a st = 1.5. 
Relevant parameters for the standard afterglow component are £j = 10 53 ' 5 erg, 
e e f = 10~ 2 , eg/ = 10 -3 , n = 1 cm 4 , and p = 2.4. For comparison, we also show 
the assumed synchrotron spectrum and the resulting SSC spectrum. Thick 
curves represent cases where the EBL attenuation is taken into account, while 
thin ones do not. Note that the attenuation by pair creation in the source is 
considered. The Feraii/LAT and CTA sensitivities (with the duty factor of 
30 %) are also overlayed (CTA Consortium 2010). 



low-energy cutoff due to synchr otron self-absorption i s given 
from the blackbody limit (e.g., IShen & Zhang 20 09). But, 
note that those cutoff energies are not relevant for our results. 
Jet opening angles of both the jets are set to 9j = 0.2 and their 
bulk Lorentz factors are assumed to be larger than 1 / 9j. 

In this section, we discuss the results on high-energy af- 
terglow emission, that is, EIC and SSC components from the 
early jet. One should keep in mind that SSC emission from 
the late jet, which was discussed in the previous section, may 
also exist. In the late internal shock and magnetic dissipation 
scenarios, one could expect ~ 1 — 10 GeV gamma rays via the 
SSC mechanism, which are potentially important for Fermi 
if nearby and/or energetic GRBs occur. But its predictions 
are highly model dependent, and very high energy emission 
is not expected when r, and/or r em are small enough (e.g., 
in the dissipative photosphere scenario), so that it will not be 
shown here. On the other hand, the EIC emission considered 
here is independent of various late internal dissipation scenar- 
ios. Even if the SSC emission from the late jet exists, this 
EIC and/or SSC components from the early jet will typically 
be dominant at very high energies. Therefore, our results on 
the EIC emission provide the most conservative high-energy 
predictions of the late internal dissipation model. 

The resulting spectra for our typical parameter sets are 
shown in Figures 3 and 4. As expected in the previous sec- 
tion, the EIC peak is located at E^ IC ~ 10-100 GeV. In our 
cases, the EIC peak energy is comparable to the SSC peak en- 
ergy at T ~ T a , which can be understood from E^ IC /E^ SC ~ 
E b /2E C . When the EIC emission is dominant, its spectrum is 
roughly expressed by Equations (21). (When the SSC emis- 
sion is dominant, its spectrum is roughly expressed by Equa- 
tion (29).) As expected before, the KN suppression becomes 
important above ~ 1 — 10 TeV but it is difficult to be observed 
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FIG. 4. — Spectra of early and late jets in the late dissipation model, con- 
sidered in this work. Syn and SSC come from synchrotron and SSC emis- 
sion by relativistic electrons accelerated at the external shock of the early 
jet. LP represents the assumed seed photon spectrum from the late jet, which 
is responsible for shallow-decay X-ray emission, and EIC is the EIC emis- 
sion by Compton scatterings of X-ray photons by electrons accelerated at the 
external shock. The observation time is set to T = 10 3 ' 6 s and the source 
redshift is taken as z = 0.3. Here, relevant parameters for the late jet are 
L*p| r = 10 48 erg s -1 , T a = 10 3 s, E' b = 0.1 keV, a s = 0.2, and a s t = 1.5. 
Parameters for the standard afterglow component are the same as those used 
in the caption of Figure 3. Thick curves represent cases where the EBL atten- 
uation is taken into account, while thin ones do not. The attenuation by pair 
creation in the source is also considered. 
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FIG. 5. — Electron cooling timescales by the EIC, SSC, and synchrotron 
processes at the external shock radius of r — 10 17 ' 5 cm are shown. For com- 
parison, the dynamical timescale is also shown. One can see that the corre- 
sponding 7 e , c ~ 10 4 5 . Source parameters are the same as those used in the 
caption of Figure 4. 



due to the EBL attenuation. 

For these parameter sets of £™p X /e e f£k ~ 1 and p ~ 2.4, 
we see that the EIC flux becomes larger than the SSC flux 
at T ~ T a . This can be understood by comparing electron 
cooling timescales. An example of three timescales is shown 
in Figure 5, where we can see that % iC ~ 10 4-5 at r ~ 10 17 ' 5 cm 
(or T - 10 3 s). When the EIC and SSC cooling times can be 
estimated in the Thomson limit, we obtain 
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FIG. 6. — Gamma-ray light curves of EIC emission at 1 GeV in the late 
internal dissipation model for GRB afterglows, caused by Compton scatter- 
ings of X-ray photons by electrons accelerated at the external shock. For 
comparison, light curves of synchrotron and SSC emission are also shown. 
The parameter set is the same as that used in the caption of Figure 3. The 
Fermi/LAT sensitivity is overlaid. Note that the attenuation by pair creation 
both inside and outside the source is taken into account. 
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FIG. 7. — Same as Figure 6, but at 100 GeV. The CTA sensitivity (with the 
duty factor of 30 %) is overlaid instead of the Fermi one (CTA Consortium 
2010). 



Then, the ratio of the EIC flux to the SSC flux is roughly 
estimated as ~ x?sscAeic- Note that, for sufficiently large j e , c 
and/or E h IE c , the results are affected by the KN effect. 

In the case shown in Figure 5, the synchrotron cooling is 
dominant. If the EIC cooling is more important than the SSC 
cooling and the synchrotron cooling, the afterglow emission 
from the early jet is affected by the EIC cooling. However, 
this occurs only when the late prompt emission from the late 
jet is bright enough. The associated afterglow emission from 
the early jet, produced by electrons with ~ j e>c , is typically 
masked by the emission from the late jet, so that it seems diffi- 
cult to observe the EIC influence at the optical or X-ray band. 

The resulting light curves are shown in Figures 6, 7 and 
8. The SSC flux evolves as EF SSC oc T~ p l 2 at E > £| sc . On 
the other hand, the EIC flux has shallower light curves, but 
its time evolution is different from that in the X-ray band 
(see Figure 8). In this sense, the EIC emission in the late 
internal dissipation model can be distinguished from the pre- 
dictions of other models, such as SSC emission from the 
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FIG. 8. — Light curves of early and late jets in the late dissipation model at 
various energy bands. Syn and SSC come from synchrotron and SSC emis- 
sion by relativistic electrons accelerated at the external shock of the early jet. 
LP represents the assumed seed photon emission from the late jet, which is 
responsible for shallow-decay X-ray emission, and EIC is the EIC emission 
by Compton scatterings of late prompt photons by electrons accelerated at 
the external shock. The parameter set is the same as that used in the caption 
of Figure 4. Note that the attenuation by pair creation both inside and outside 
the source is taken into account. 



late jet or SSC afterglow emission in modified FS mod- 
els. The time evolution of the EIC emission is understood 
from F E ic = t^\ c {^ e ^ c )/t~^ e ^ c ). If electrons with -f ex are 
in the Thomson regime (which is not always true), we ex- 
pect Ye IC oc L lp /R 2 T 2 B 2 oc r _QLP+I . On the other hand, 
the synchrotron luminosity in the slow cooling case obeys 
^ssc ^ T~ p l 2 from Equation (30). Then, we roughly expect 
^eic^eic ~ xY mc {E c F c ) oc r- QLP+I -''/ 2 , which declines more 
rapidly than the shallow decay emission. For example, for 
p ~ 2.4 and a LP ~ 0.2, we have EF E i C oc T~ 0A . When the KN 
effect plays a role, the temporal index is somewhat steeper, 
which seems consistent with the numerical results. The break 
time of the shallow-decay emission is T a ~ 10 3 s, but the EIC 
flux does not decline for a while even after T„. This is be- 
cause seed photons interacting with FS electrons come from 
backward (9 ~ 0), leading to suppressed EIC emission toward 
the observer, while photons passing through the FS with an- 
gles of 9 y 0, which significantly contribute to the EIC flux, 
are delayed compared to non-scattered photons from 9 ~ 0. 
Note that the timescale of this delay is of order of R/T 2 c ~ T, 
which is understood from the fact that EIC emission induced 
by an impulsive seed photon emission lasts until we observe 
photons entering the FS region with 9 ~ 1/r. 

From Figures 3 and 7, for our typical parameter sets, the 
EIC emission is expected at energies larger than 10 GeV. 
As we can see, the EBL attenuation is moderate for nearby 
GRBs, though it becomes crucial for distant bursts (see be- 
low). At such very high energies, observations by Cherenkov 
telescopes such as MAGIC, VERITAS, HAWC, and CTA 
are more promising. Although no clear detections have 



been obtained so far dAbdo et al 
Aharonia net al. 20091; 



2007; 



Albert et al. 2001; 



Aleksic et al. 2010). future observa 



tions with HAWC and CTA would improve the chances for 
this, and either detections or non-detections are important to 
test the model. Detections by Fermi are limited at late times, 
but they are being made in the earlier afterglow phase. How- 
ever, note that the synchrotron or SSC emission is more im- 
portant than the EIC emission at the earlier phase (especially 
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FIG. 9. — Gamma-ray spectra of EIC and SSC emission in the late internal 
dissipation emission model at T = 10 2 6 s (thin) and T = 10 3 6 s (thick), but 
the source redshift is taken as z = 1 ■ Source parameters are the same as that 
used in the caption of Figure 3. The EBL attenuation is included in EIC and 
SSC, but not included in EIC (intrinsic) and SSC (intrinsic). One can see that 
it is crucial for detections by Cherenkov telescopes. 
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FIG. 10. — Gamma-ray light curves of EIC emission at 1 GeV in the late 
internal dissipation model for GRB afterglows. For comparison, light curves 
of synchrotron and SSC emission are also shown. The parameter set is the 
same as that used for Figure 9. Here, the attenuation by pair creation both 
inside and outside the source is taken into account. 



just after the prompt emission), which can be expected from 
Equation (40). Also, the SSC emission from the late jet, 
which has been discussed in the previous section, can be rel- 
evant in the GeV range (e.g., in the magnetic dissipation sce- 
nario), while the EIC emission and/or SSC emission from the 
early jet will be dominant in the 100 GeV range. 

We have considered in our calculations the EBL attenu- 
ation, using the low-IR model developed by Kneiske et al. 
(2004). Detecting gamma rays at very high energies above 
100 GeV is prevented by this EBL attenuation. Even for a 
burst at z = 0.3, we have seen that the EBL attenuation largely 
degrades the resulting fluxes at > 300 GeV. For a burst at 
higher redshifts, the situation becomes worse. In Figures 9, 
10 and 11, results for z = 1 are shown, where the EBL attenu- 
ation becomes crucial at > 100 GeV. The EIC and SSC peak 
energies are higher than the cutoff by the EBL attenuation, 
so that both of the EIC and SSC fluxes are largely degraded. 
However, detections around ^10 GeV appear still promising 
at earlier times, even though they are difficult at late times. 
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FIG. 11 .— Same as Figure 10 but at 100 GeV. 



Note that gamma rays absorbed by the EBL must produce 
energetic pairs, which lead to IC emission by scatterings with 
the EBL photons. In our calculations, this pair echo emis- 
sion is not included since it is beyond the scope of this work, 
although it could affect the observed afterglow emission if 
the intergalactic magnetic field in voids i s weak enough (e.g., 
iRazzaque et al. 2004tlMurase et al. 20091) . 

4.1. Discussion on Parameter-Dependence 

We have demonstrated that the EIC emission dominates 
over the SSC emission in the late internal dissipation model 
for shallow-decay emission. Importantly, predictions of the 
EIC emission are straightforward, once X-ray and optical af- 
terglows are well observed. The parameters necessary for cal- 
culations of high-energy emission are determined via fitting 
with the two-component (early and late) jet model as done in 
Ghisellini et al. (2009). The parameter dependence of the 
relative importance of the EIC emission to the SSC emission 
is seen from Equation (1401) . The most important quantity is 
£^p/e e f£k (which is expected by setting p ~ 2), which is seen 
by comparison between Figures 7 and 8. For our typical cases, 
the EIC emission is dominant at late times (T > T a ), but can 
be less important for smaller values. In fact, there is large di- 
versity among observed X-ray and optical afterglows so that 
it would be natural to expect that high-energy afterglows also 
exhibit a high diversity, depending on £ L s p / e e /£k. 

In our calculations, we have assumed egf = 10~ 3 , but the 
EIC and SSC peaks are rather sensitive to e^/ (see Equa- 

so that 



£ 3/2 
£ Bf 



tions (17) and (28)). We see that ^eic/^ssc 
the EIC peak is more likely to be higher than the SSC peak 
for larger eg/. This implies that the EIC component is more 
frequently dominant over the SSC one at high energies. Note 
that the EIC peak energy can be around 1 — 10 GeV rather than 
0.1-1 TeV when e Bf ~ 10" 2 . 

Another potentially relevant parameter is E b . For typical 
values used in this work, the results on the EIC emission are 
not so sensitive to this quantity, up to a modest factor (see Fig- 
ure 12). But this may not be the case if the EIC cooling occurs 
in the KN regime. If E b is so large that the EIC cooling occurs 
in the KN regime while the SSC cooling does in the Thom- 
son regime, the EIC emission would be more suppressed. So 
far, we have assumed that E b does not depend on time. Al- 
though this may not be true, it is difficult to determine its time 




log(E y [GeV]) 

FIG. 12. — EIC spectra calculated with different values of E b . Intrinsic EIC 
spectra (where the EBL attenuation is not taken into account) are shown at 
T = I0 2 6 s (thin) and T = 10 3 ' 6 s (thick). The source redshift is set to z = 0.3. 
For the solid curves, the used parameter set is the same as that in the caption 
of Figure 3. The dotted curves are for E' b = 10 2 5 eV and the dashed curves 



are for E' b = 0.1 keV (L* p /L* p | r 



1/2 



while the other parameters are the 
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FIG. 13. — EIC light curves (thick) at 100 GeV, calculated with different 
assumptions on E h and IYp. For comparison, SSC light curves (thin) are 
also shown. The source redshift is set to z = 0.3. For the solid curves, the 
used parameter set is the same as that in the caption of Figure 3. The dashed 

curves are for E' b = 0.1 keV (L* p /L* p | To /)' /2 rather than E' b =const. The 
dotted curves are for a s t = 2.8 rather than a s t = 1-5. The dot-dashed curves 
are for the case with the evolving E' and a st = 2.8. The attenuation by pair 
creation both inside and outside the source is taken into account. 



evolution from observations. To see the influence of this un- 
certainty on results, we also calculate the EIC emission with 

the break energy of E b {T) oc L\£{T). However, as seen from 
Figures 12 and 13, the results are hardly changed, because the 
EIC emission mainly occurs in the Thomson regime. At later 

times, the EIC flux with E b (T) oc l\^(T) is a bit larger than 
that with E b =const. This is because lower E b at late times 
can compensate the KN effect due to increase of 7 e c . 

The low-energy photon index of the late prompt emission, 
Pi, is observationally uncertain. But this becomes crucial for 
the EIC spectrum at relatively low energies of E < E£ 1C , so 
that we expect that our results are not affected by this. In 
addition, the EIC flux is larger if /?/ > 1 . 

We also show, for comparison, the resulting EIC and SSC 
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light curves for a st = 2.8 in Figure 13, where one can see that 
the EIC light curve declines more rapidly after T a . For the jet 
opening angle and bulk Lorentz factor of the late jet, we have 
assumed Tlp > 1 / $lp- This may not be the case, as discussed 
in Ghisellini et al. (2007). If the late jet is decelerated with 
time, we expect that the observed light curve has the break 
when Tlp becomes ~ 1/#lp- This break may be the origin 
of t a , although it is not clear why the late jet is decelerated 
continuously. In this case, only a fraction of seed photons can 
interact with FS electrons after T a . But this just corresponds 
to a change of a, which is already taken into account in ob- 
servable parameters. 

The jet opening angles and axes of the two jets are also as- 
sumed to be the same. However, we can still expect the EIC 
emission even if the jets are a bit misaligned. If either edge of 
the early jet is on the line of sight, photons from the late jet 
still come to the observer through the early jet (independently 
of the prompt emission mechanism), but the resulting EIC flux 
is reduced by a factor of two at most. Note that the important 
assumption used in this work is (r,-/ r)T6up "C 1, which is typ- 
ically valid in our model. When this condition does not hold, 
more detailed calculations are required. 

4.2. Specific Cases 

Bursts with a sudden decline in their X-ray afterglows 
are of particular interest. For example, GRB 070110 
has a steep decline o f a st = 9 after the plateau of an = 
0.09 dTroia et al. 20071) . In Figures 14 and 15, we show the 
specific case of GRB afterglows with such a plateau (with 
Qffl = and a st = 10), to see the EIC emission induced by the 
plateau X-ray emission. It is obvious that the EIC spectrum is 
similar to that shown in Figure 3, since a similar seed photon 
spectrum is assumed. On the other hand, the EIC light curve 
is different from those shown in Figure 4, reflecting differ- 
ent X-ray light curves. Before T ~ T a , the EIC light curve is 
steeper than the X-ray one, as discussed before. However, this 
is not the case after ~ T a . This is because the EIC emission 
is similar to high-latitude emission, so that the EIC emission 
does not show a sudden decline even though the seed photon 
emission ends abruptly. This was the behavior seen for an 
impulsive seed p hoton emission, as demonstrated in the prior 
emission model dMurase et al. 20101) . Detecting such a sig- 
nature of high-latitude emission associated with the sudden 
decline after the plateau would be useful as evidence of the 
late internal dissipation model. 

Another case of interest is that of the GRBs that were ob- 
served by Fermi which may be represented by the late dissi- 
pation models. The high-energy emission detected by Fermi 
may originate from the external shock. For example, Kumar 
& Barniol Duran (2010) argued that long-lasting GeV emis- 
sion comes from electrons accelerated at the FS caused by an 
adiabatic relativistic blast wave expanding into a low density 
ISM (n — 10" 4 cm" 3 ), with a low magnetic field (e B f ~ 10" 4 ). 
On the other hand, Ghisellini et al. (2010) argued that long- 
lasting GeV emission may be explained by a radiative rel- 
ativistic blast wave, with e e f ~ 1 and p ~ 2. Although the 
origin of GeV emis sion especially at the very early stage 
is still under debate dHe et al. 20101; iLiu & Wang 20111) , the 
late-time GeV emission is likely to be regarded as afterglows. 

Just for demonstrative purposes, we also calculate the EIC 
emission for a burst like GRB 090902B. Unfortunately, there 
have been no bursts that have canonical early afterglow light 
curves simultaneously observed by Fermi and Swift, so that 
we just show the result for parameters provided in Cenko et al. 



SSC 




FIG. 14. — Gamma-ray spectra of EIC emission from the GRB afterglow 
with the plateau X-ray emission at T = 10 2 6 s (thin) and T = 10 3 6 s (thick). 
For comparison, synchrotron and SSC emission from the standard afterglow 
component are also shown. The source redshift is set to z = 0.3. Relevant pa- 
rameters for the late jet are z£ p | T , = 10 48 5 ergs" 1 , T' a = 10 3 s, £'* = 0.1 keV, 
Ofl = 0, and a st = 10. Relevant parameters for the standard afterglow compo- 
nent are the same as those in the caption of Figure 3. The attenuation by pair 
creation both inside and outside the source is taken into account. 
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FIG . 15. — Light curves of early and late jets in the late dissipation model at 
various energy bands. Syn and SSC come from synchrotron and SSC emis- 
sion by relativistic electrons accelerated at the external shock of the early jet. 
LP represents the assumed seed plateau emission from the late jet, and EIC is 
the EIC emission by Compton scatterings of late prompt photons by electrons 
accelerated at the external shock. The parameter set is the same as that used 
for Figure 14. We can see that the EIC light curve is much shallower than 
that of late prompt emission after T a . The attenuation by pair creation both 
inside and outside the source is taken into account. 



(2010) in Figure 16 (see also Pandey et al. 2010) (but the red- 
shift is set to z = 1 here). Those parameters are indicated from 
late-time observations but not exact ones for explaining long- 
lasting GeV emission, since implied j e>c is smaller than those 
used in Kumar & Barniol Duran (2010). In fact, Liu & Wang 
(2010) showed that the GeV emission may rather be explained 
by an additional jet component (although the two-component 
jet model used there is different from that considered here). 
Here, we are not pursuing possibilities to explain the long- 
lasting GeV emission with the late internal dissipation model, 
so these parameters are suitable enough for the present pur- 
pose. As can be seen, for our parameters on late prompt 
emission, the EIC component is dominant over the SSC one 
at late times. The EIC component is especially important at 
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FIG. 16. — Gamma-ray spectra of EIC emission for a burst with GRB 
090902B-like afterglow parameters, at T = 10 2 6 s (thin) and T = 10 36 s 
(thick). For comparison, synchrotron and SSC emission from the standard 
afterglow component are also shown. The source redshift is set to z = 1 ■ Rel- 
evant parameters for the late jet are L^p\t' = 1-02 X 10 50 erg s _1 , T' a = 10 3 s, 
E' =0.1 keV, «a = 0, and a st = 2. Relevant parameters for the stan- 
dard afterglow component are Sk = 6.8 X 10 53 erg, e e f = 0.15, esf = 0.058, 
n = 5.8 X 10 -4 cm -3 , and p = 2.22. The attenuation by pair creation both 
inside and outside the source is taken into account. 



high energies of ~ 10-100 GeV. At lower energies, the syn- 
chrotron component dominates over the others, although the 

curve shown is fairly optimistic (since j e ,M = \J is used 

assuming that the upstream magnetic field is the downstream 
one and r\ ~ 1). As is demonstrated in Figure 16, the EIC 
emission induced by late prompt emission can typically be 
important in very high energies only at relatively late times, 
so that our results on high-energy afterglow emission are com- 
patible with the Fermi observations. On the other hand, high- 
energy late prompt emission could potentially be relevant in 
the GeV range. One may expect that it shows the shallow- 
decay behavior when the shallow-decay emission comes from 
the late internal dissipation. However, the observational sit- 
uation is currently unclear since simultaneous detections by 
Fermi and Swift are required. Possibly, for bursts detected 
by FermifLP^, T a is large enough and it may become impor- 
tant only at late times, or T a is small enough but it may be 
masked. Or, declining high-energy emission could happen, 
if the steep-decay emission comes from the late internal dis- 
sipation. Also, ferm//LAT bursts tend to be most energetic 
ones, and it has not been settled whethe r accelerated electrons 
are in the fast or slow coolin g regime (Ghis elhni et al. 20101: 
iKumar & Barniol Duran 2 010). Future simultaneous detec- 
tions of high-energy gamma rays from GRBs with canonical 
afterglow light curves are anticipated. 

5. SUMMARY AND DISCUSSION 

In this paper, we have studied the possibility that the 
shallow-decay or plateau emission originates from late inter- 
nal dissipation in the late jet driven by the long-lasting central 
engine (e.g., mass fall back onto a black hole or rotational 
energy loss of fast rotating magnetars). We have discussed 
various theoretical scenarios of the emission mechanism, late 
internal shock, magnetic dissipation, and photospheric sce- 
narios. There are few clues to the origin of the late prompt 
emission, and all the three scenarios seem compatible with 
observations at present. 



We have also investigated the associated high-energy emis- 
sion in the late internal dissipation model and discussed two 
possibilities: high-energy late prompt emission and high- 
energy afterglow emission. The former comes from internal 
dissipation in the late jet, and the predictions depend on the 
specific scenarios. For example, in the photospheric scenario, 
high-energy emission may be produced by the IC process but 
> GeV emission is not expected in the one-zone case due to 
the large pair-creation opacity around the photosphere. On 
the other hand, the late internal shock and magnetic dissipa- 
tion scenarios may lead to ~ 1 - 10 GeV gamma rays by the 
SSC mechanism, which are important for Fermi. As demon- 
strated in this work, detections by Fermi and possibly CTA are 
expected for nearby and/or energetic GRBs, which would be 
useful for revealing the mechanism of late prompt emission. 

The latter possibility includes the SSC and EIC emis- 
sions produced by electrons accelerated at the external shock, 
which will be especially relevant in the very high energy 
range. Especially, the EIC emission, which is high-energy 
afterglow emission induced by late prompt photons, is not so 
sensitive to details of models and should be useful as a test of 
the existence of late internal dissipation during the shallow- 
decay phase. In this work, we have investigated the EIC emis- 
sion both analytically and numerically and demonstrated that 
the EIC flux may become larger than the SSC flux around 
the end time of the shallow-decay phase. The EIC peak is 
typically expected at ~ 1-100 GeV, and the EIC emission 
typically has a steeper light curve than in the X-ray one, but 
a shallower one when the X-ray light curve shows a sudden 
decline. Hence, it would be possible to distinguish it from the 
other possibilities such as SSC components from the early and 
late jets. We also expect that it is easier for the synchrotron 
and SSC components to dominate at very earlier times. 

Although the detectability depends on the parameters 
and on the EBL, ground-based gamma-ray observatories, 
such as MAGIC, VERITAS, HESS, CTA and HAWC, 
would be important tools in the search for such signals. 
Very high energy gamm a rays from GRBs have not been 
firmly observed so far dAbdo et al. 20071: lAlbert et al. 20071: 
Aharo riian et al. 2009t lAleksic et al. 2010h and the event rate 
of nearby bursts is not large. (For example, the rate of GRBs 
occurring at within z ~ 0.3 is es timated as ~ a few events 
per year (e.g. jLiang et al. 2007alU Nevertheless, once suffi- 
ciently fast follow-up observations are successful for nearby 
events, Cherenkov telescopes with a low-energy threshold 
(~ 30 GeV) may allow us to have good photon statistics 
thanks to their high sensitivities. Theoretical predictions of 
the EIC emission are testable once parameters are specified 
from observations, and the strategy for testing the model is as 
follows. First, one determines the relevant standard afterglow 
parameters for the early flow. When afterglows are well ob- 
served at optical (and/or X-ray) bands, the parameters such as 
£k> P< e Bf and e e f are determined in the context of the stan- 
dard external FS theory. At the same time, parameters on the 
late prompt emission, such as L* P and E h and T a , can also be 
determined from observations at X-ray (and/or optical) bands. 
With those parameters, both the EIC and SSC emissions are 
calculated and can be compared to high-energy observations. 
Even non-detections would provide useful constraints on the 
models, especially for GRBs with a strong plateau or shallow- 
decay emission. 

It is also important to keep in mind that GRB after- 
glows seem to be fairly diverse (e.g., iGhisellini et al. 20091: 
iLiang et al. 20091) . Although GRBs with a shallow-decay 
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emission may be explained by a late jet from the long-lasting 
central engine, some GRBs do not show the shallow decay 
and can be explained by the standard afterglow model, where 
the synchrotron or SSC emission is expected to be dominant. 
For this reason, multi-wavelength observations from radio to 
gamma rays are important for comprehensive studies of GRB 
afterglows. 

6. ACKNOWLEDGMENTS 

K.M. acknowledges financial support by a Grant-in-Aid 
from JSPS, from CCAPP and from PSU. K.T. and P.M. 



acknowledge partial support from NASA NNX08AL40G, 
NASA NNXAT72G, NSF PHY-0757155 and U.R.A. 10-S- 
017. This research was also supported by Grant-in-Aid from 
the Ministry of Education, Culture, Sports, Science and Tech- 
nology MEXT of Japan, no. 19047004 and 21740184 (R.Y.). 
The numerical calculations were carried out on Altix3700 
BX2 at YITP in Kyoto University. P.M. acknowledges the 
hospitality of the Institute for Advanced Study, Princeton, 
during part of this project. 



REFERENCES 



Abdo, A. A., et al. 2007, ApJ, 666, 361 
Abdo, A. A., et al. 2009a, Science, 323, 1688 
Abdo, A. A., et al. 2009b, ApJ, 706, L138 
Ackermann, M., et al. 2010, ApJ, 712, 558 
Aharonian, F. A., & Atoyan, A. M. 1981, Ap&SS, 79, 321 
Aharonian, F. A., et al. 2009, A&A, 495, 505 
Albert, J. et al. 2007, ApJ, 667, 358 
Aleksic, J., et al. 2010, A&A, 517, A5 
Baring, M. G. 2006, ApJ, 650, 1004 
Beloborodov, A. M. 2005, ApJ, 618, L13 
Beloborodov, A. M. 2010, MNRAS, 407, 1033 
Beloborodov, A. M., et al. 2011, MNRAS, 410, 2422 
Blandford, R. D., & McKee, C. F. 1976, Phys. Fluids, 19, 1130 
Blumenthal, G. R., & Gould, R. J. 1970, Rev. Mod. Phys., 42, 237 
Bosnjak, Z., Daigne, F, & Dubus, G. 2009, A&A, 498, 677 
Bottcher, M., & Dermer, C. D. 1998, ApJ, 499, L131 
Cenko, S. B., et al. 2010, arXiv: 1004.2900 
Chincarini, G., et al. 2010, MNRAS, 406, 2113 
CTA Consortium 2010, arXiv: 1008.3703 
Dai, Z. G., & Lu, T. 1998, A&A, 333, L87 
Daigne, F, Bosnjak, F, & Dubus, G. 2011, A&A, 526, A110 
Derishev, E. V., Kocharovsky, V. V., & Kocharovsky, VI. V. 2001, A&A, 372, 
1071 

Dermer, C. D. 2007, ApJ, 664, 384 

Dermer, C. D., Chiang, J., & Mitman, K. 2000, ApJ, 537, 785 
Eichler, D., & Granot, J. 2006, ApJ 641, L5 
Eichler, D., & Waxman, E. 2005, ApJ, 627, 861 
Falcone, A. D. et al. 2007, ApJ, 671, 1921 

Fan, Y. Z., & Piran, T. 2008, Frontiers of Physics in China, 3, 306 

Fan, Y. Z., Piran, T., Narayan, R., & Wei, D. M. 2008, MNRAS, 384, 1483 

Fan, Y. Z., & Wei, D. M. 2005, MNRAS, 364, L42 

Genet, F, Daigne, F, & Mochkovitch, R. 2007, MNRAS, 381, 732 

Ghisellini, G., et al. 2007, ApJ, 658, L75 

Ghisellini, G., et al. 2009, MNRAS, 393, 253 

Ghisellini, G., et al. 2010, MNRAS, 403, 926 

Giannios, D. 2010, MNRAS, 408, L46 

Giannios, D., et al. 2009, MNRAS, 395, L29 

Gou, L. J., & Meszaros, P. 2007, ApJ, 668, 392 

Granot, J., Cohen-Tanugi, J., do Couto e Silva, E. 2008, ApJ, 677, 92 

Granot, J., Piran, T., & Sari, R. 1999, ApJ, 513, 679 

Guetta, D., & Granot, J. 2003, MNRAS, 340, 115 

Gupta, N., & Zhang, B. 2008, MNRAS, 384, LI 1 

He, H. N., et al. 2010, arXiv: 1009. 1432 

Ioka, K. 2010, Prog. Theo. Phys., 124, 667 

Ioka, K., Kobayashi, S., & Zhang, B. 2005, ApJ, 631, 429 

Ioka, K., Toma, K., Yamazaki, R., & Nakamura, T. 2006, A&A, 458, 7 

Ioka, K., et al. 2007, ApJ, 670, L77 

Jin, Z. P., Fan, Y. Z., & Wei, D. M., 2010, ApJ, 724, 861 

Kneiske, T. M., et al. 2004, A&A, 413, 8075 

Kumar, P., & Barniol Duran, R. 2010, MNRAS, 409, 226 

Kumar, P., Narayan, R., & Johnson, J. L. 2008a, MNRAS, 388, 1729 

Kumar, P., Narayan, R., & Johnson, J. L. 2008b, Science, 321, 376 

Kumar, P., & Panaitescu, A. 2000, ApJ, 541, L51 

Lazar, A., Nakar, E., & Piran, T. 2009, ApJ, 695, L10 

Liang, E. W., Zhang, B., Virgili, F, & Dai, Z. G. 2007a, ApJ, 662, 1111 

Liang, E. W., Zhang, B. B., & Zhang, B. 2007b, ApJ, 670, 565 

Liang, E. W., et al. 2009, ApJ, 707, 328 



Lithwick, Y, & Sari, R. 2001, ApJ, 555, 540 
Liu, R. Y. & Wang, X. Y. 201 1, ApJ, 730, 1 

Lyutikov, M. 2006, New J. Phys., 8, 1 19 

McKinney, J. C, & Uzdensky, D. A. 2010, arXiv:1011.1904 
Meszaros, P. 2006, Rep. Prog. Phys., 69, 2259 
Meszaros, P., & Rees, M. J. 1997, ApJ, 476, 232 
Meszaros, P., & Rees, M. J. 2000, ApJ, 530, 292 

Meszaros, P., & Ramirez-Ruiz, E., Rees, M. J., & Zhang, B. 2002, ApJ, 578, 

812 

Metzger, B. D., et al. 2010, arXiv: 1012.0001 

Murase, K. 2007, Phys. Rev. D, 76, 123001 

Murase, K, & Ioka, K. 2008, ApJ, 676, 1123 

Murase, K, & Nagataki, S. 2006, Phys. Rev. Lett., 97, 051101 

Murase, K, Zhang, B., Takahashi, & Nagataki, S. 2009, MNRAS, 396, 1825 

Murase, K, et al. 2010, MNRAS, 402, L54 

Nakar, E., Ando, S., & Sari, R. 2009, ApJ, 703, 675 

Nousek, J. A., et al. 2006, ApJ, 642, 389 

O'Brien, P.-T., et al. 2006, ApJ, 647, 1213 

Panaitescu, A. 2008a, MNRAS, 383, 1143 

Panaitescu, A. 2008b, MNRAS, 385, 1628 

Panaitescu, A., & Kumar, P. 2004, MNRAS, 353, 511 

Panaitescu, A., et al. 2006, MNRAS, 366, 1357 

Pandey, S. B., et al. 2010, ApJ, 714, 799 

Pe'er, A., Meszaros, P., & Rees, M. J. 2006, ApJ, 642, 995 

Pe'er, A., & Waxman, E. 2005, ApJ, 633, 1018 

Razzaque, S., Meszaros, P., & Zhang, B. 2004, ApJ, 613, 1072 

Rees, M. J., & Meszaros, P. 1994, ApJ, 430, L93 

Rees, M. J., & Meszaros, P. 1998, ApJ, 496, LI 

Rees, M. J., & Meszaros, P. 2005, ApJ, 628, 847 

Rhoads, J. E. 1999, ApJ, 525, 737 

Sari, R., & Esin, A. A., 2001, ApJ, 548, 787 

Sari, R., Piran, T., & Halpern, J. 1999, ApJ, 519, L17 

Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 

Shen, R. F, & Zhang, B. 2009, MNRAS, 398, 1936 

Thompson, C. 1994, MNRAS, 270, 480 

Thompson, T. A., Chang, P., & Quataert, E. 2004, ApJ, 611, 380 

Toma, K, Ioka, K„ Yamazaki, R., & Nakamura, T. 2006, ApJ, 640, L139 

Toma, K. Wu, X. F, & Meszaros, P. 2009, ApJ, 707, 1404 

Troja, E., et al. 2007, ApJ, 665, 599 

Uhm, Z. L., & Beloborodov, A. M. 2007, ApJ, 665, L93 

Wang, X. Y, Li, Z., & Meszaros, P. 2006, ApJ, 641, L89 

Wang, X. Y, et al. 2010, ApJ, 712, 1232 

Willingale, R., et al. 2007, ApJ, 662, 1093 

Woods, E., & Loeb, A. 1999, ApJ, 523, 187 

Yamazaki, R. 2009, ApJ, 690, LI 18 

Yamazaki, R., Toma, K, Ioka, K., & Nakamura, T. 2006, MNRAS, 369, 311 

Yu, Y. W., Cheng, K. S., & Cao, X. F. 2010, ApJ, 715, 477 

Zhang, B. 2007, Chin. J. Astron. Astrophys., 7, 1 

Zhang, B., & Meszaros, P. 2001a, ApJ, 552, L35 

Zhang, B., & Meszaros, P. 2001b, ApJ, 559, 110 

Zhang, B., & Yan, H. 201 1, ApJ, 726, 90 

Zhang, B., et al. 2006, ApJ, 642, 354 

Zhang, B., et al. 2007, ApJ, 655, 989 

Zhang, B. B., et al. 2009, ApJ, 690, L10 



High-Energy GRB Afterglows in the Early and Late Jet Model 



15 



APPENDIX 
FORMULAS OF EIC EMISSION 



Here, we de rive formulas of EIC emission in GRB aft erglows. The observed flux from the shell expanding toward us relativis- 
tically is (e.g., iGranotetal. 19991; I Woods & Loeb 19991) 



F(7) = i±! / d« j </cosf 



drr 1 



.h 



r 2 ( i-/3 cos ey 



(Al) 



where j £ is the comoving emissivity and e is the seed photon energy in the comoving frame. Hereafter, we also use E = 
(1 +z)ET(\— j3cos9) and T = (l+z)(f — rcos8/c). Especially, the comoving EIC emissivity is written as (e.g jToma et al. 20091) 
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Here, y : 



£m e c" 



■, $ = — — ~ , and scattering angles 9 and 9 are measured in the central engine frame and the comoving 
frame, respectively. The range of y is t— — J ln _ e -, < y < 1 ■ Note that Equation (A2) is easily obtained from 



dNl 
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(A4) 



where b § = 2(1- cos 9)^f e e/m e c 2 dAharonian & Atovan 198ltlFan & Piran 20081) . 

First, we shall derive the formula for an impulsive seed photon spectrum. We also assume that seed photons come from 
r,- r. In the case of instantaneous emission (at to) from an infinitely thin shell (at Ro), by using the replacement of j £ — > 
j E 6(f-T ) f dyn 5(r-R Q ) A, we obtain 

e 



F Elc (T)=-<T T (l-cos9) / d% 



dn e A 

d~f e K 



dy Fseed|r ° , (1-0 

(i+r 2 e 2 ) 2 
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where fdyn = A/c is the comoving dynamical timescale and A = Rq/kT is the comoving shell thickness. Here, 



2 (D = 2 



1 [To 
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In the case of a broken power-law seed spectrum, we can write F se ed|r = Fs eed \ t0 G(e), where ^ S e e dl7o = ^jrjk^ (1+z) which is 
smeared over the dynamical timescale of the shell. Note that AT is the duration of impulsive seed photon emission in the oberver 
frame. Equation (A5) is essentially the same as Equation (5) used in Murase et al. (2010) 1 . 

Next, we shall derive the formula for continuous seed photon emission. This is obtained by the similar procedure. Performing 
the replacement of j e —> j £ 5(r-R(T)) A(T) leads to 



Feic(T)= -err J dr (I -cost 
where A = r/nT and 9 = 9(r) is given by 
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and 9 = 6(f) is obtained via the Lorentz transformation. When T9^$> 1, we obtain Equation 

DISTRIBUTION OF NONTHERMAL ELECTRONS 

In order to calculate both the EIC and SSC emission, we use the following electron distribution for j e > j ejn , which would 
approximately mimic the distribution of relativistic electrons in the dynamical timescale, 



dn e 

d^ e 



cxmin[lX 1 0l ]7/, 



(Bl) 



where p is the spectral index of accelerated electrons and / coo i = fd yn Acooi is the effective optical depth for energy losses. In the 
slow cooling case with f coo i = f syn (where f syn is the synchrotron cooling timescale), we have dn e /d^ e oc 7~ p for j e>m < j e < j e>c 



7 There was an unimportant typo in that paper, but calculations were performed using the correct expression, dropping off k. 



16 



Murase et al. 



and dn e /d-f e oc j e p 1 for -f e > j e , c . In the fast cooling case, we set p = 1 for ^ e , c < 7 g < 7 e . m , whic h reproduces oc j e 2 if f coo i = f syn . 
The value of j eiC is determined by finding solutions of (e.g. jNakar et al. 2 009: Wa ng et al. 20101) 



C = f S yn(7.) + hscile) + tmche), 



where the IC loss timescales are evaluated from 



he ' = 



(7,-D 
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where K\c is the electron inelasticity for the IC process (which is calculated from Equation (C2)) and <j\c is the IC cross section 
which is given by the KN formula. 
The normalization is determined by 



dn e 



4tt 



aj e 3—k 



(B4) 



where k = for the ISM and k = 2 for the wind medium. (A somewhat different normalization, Af e w 47rr 2 (4r«)(r/4r) was used 
in Murase et al. 2010.) 



SSC EMISSION AND PAIR PRODUCTION 

In this work, we also calculate the SSC emission for comparison. For simplicity, we simply calculate the observed 
SS C flux from the comoving SSC power per comoving energy. The comoving SSC power per comoving energy is given 
by dBlumenthal & Gould 1970h 
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and v = 
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and w '■■ 



4e7 P 
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Note that numerically calculated SSC fluxes have convex curves, which lead to larger fluxes 



compared to analytically c alculated SSC seg ments dSari & Esin 2001b . and the Klein-Nishina effect becomes often important 
above the TeV range (e.g., Wang et al. 2010). As for the seed photon spectrum, the an alytical synchro tron spectrum is used in 
this work, which is expressed as three segments both in the fast and slow cooling cases (Sar i et al. 19981) . 

High-energy gamma rays may suffer from pair-production process with target photons in the source. We also take into account 
the resulting gamma-ray attenuation in the source. The optical depth for the pair production is expressed as 



A f f dn syn 3 7 

r 77 (£)=— / dfi(l-n) I de^-—a T (l-/3£ M ) 



de 16 



2/3cM(^ M -2) + (3-^ M )ln' ' ' CM 



CM 



(C3) 



where /3cm = yJX -4m 2 c 4 /S) and S is the Mandelstam variable. In this work, pair attenuation in the source is taken into ac- 
count by introducing the suppression factor 1/(1 +r 77 ) ([Baring 2006). More detailed discussions on the opacity effect is found 
in Granot et al. (2008). When the gamma-ray attenuation becomes crucial, one has to consider electromagnetic cascades in 
the source. However, as long as we consider the afterglow emission in the late phase as done in this work, the gamma-ray 
attenuation is not important, and we can neglect the electromagnetic cascades caused by the leptonic SSC and EIC emission. 
The situation is different wh en one consider hadronic scenarios, where hadronic cascades may be important at very high ener- 
gies ([iVir&JVaimlnJQO!). 



